Expression network analysis of bovine skin infested with Rhipicephalus australis identifies pro-inflammatory genes contributing to tick susceptibility.

The skin is the primary feeding site of ticks that infest livestock animals such as cattle. The highly specialised functions of skin at the molecular level may be a factor contributing to variation in susceptibility to tick infestation; but these remain to be well defined. The aim of this study was to investigate the bovine skin transcriptomic profiles of tick-naïve and tick-infested cattle and to uncover the gene expression networks that influence contrasting phenotypes of host resistance to ticks. RNA-Seq data was obtained from skin of Brangus cattle with high (n = 5) and low (n = 6) host resistance at 0 and 12 weeks following artificial tick challenge with Rhipicephalus australis larvae. No differentially expressed genes were detected pre-infestation between high and low resistance groups, but at 12-weeks there were 229 differentially expressed genes (DEGs; FDR < 0.05), of which 212 were the target of at least 1866 transcription factors (TFs) expressed in skin. Regulatory impact factor (RIF) analysis identified 158 significant TFs (P < 0.05) of which GRHL3, and DTX1 were also DEGs in the experiment. Gene term enrichment showed the significant TFs and DEGs were enriched in processes related to immune response and biological pathways related to host response to infectious diseases. Interferon Type 1-stimulated genes, including MX2, ISG15, MX1, OAS2 were upregulated in low host resistance steers after repeated tick challenge, suggesting dysregulated wound healing and chronic inflammatory skin processes contributing to host susceptibility to ticks. The present study provides an assessment of the bovine skin transcriptome before and after repeated tick challenge and shows that the up-regulation of pro-inflammatory genes is a prominent feature in the skin of tick-susceptible animals. In addition, the identification of transcription factors with high regulatory impact provides insights into the potentially meaningful gene-gene interactions involved in the variation of phenotypes of bovine host resistance to ticks.

Application of quantitative proteomics to discover biomarkers for tick resistance in cattle.

Introduction: Breeding for tick resistance is a sustainable alternative to control cattle ticks due to widespread resistance to acaricidal drugs and the lack of a protective vaccine. The most accurate method used to characterise the phenotype for tick resistance in field studies is the standard tick count, but this is labour-intensive and can be hazardous to the operator. Efficient genetic selection requires reliable phenotyping or biomarker(s) for accurately identifying tick-resistant cattle. Although breed-specific genes associated with tick resistance have been identified, the mechanisms behind tick resistance have not yet been fully characterised. Methods: This study applied quantitative proteomics to examine the differential abundance of serum and skin proteins using samples from naïve tick-resistant and -susceptible Brangus cattle at two-time points following tick exposure. The proteins were digested into peptides, followed by identification and quantification using sequential window acquisition of all theoretical fragment ion mass spectrometry. Results: Resistant naïve cattle had a suite of proteins associated with immune response, blood coagulation and wound healing that were significantly (adjusted P < 10- 5) more abundant compared with susceptible naïve cattle. These proteins included complement factors (C3, C4, C4a), alpha-1-acid glycoprotein (AGP), beta-2-glycoprotein-1, keratins (KRT1 & KRT3) and fibrinogens (alpha & beta). The mass spectrometry findings were validated by identifying differences in the relative abundance of selected serum proteins with ELISA. The proteins showing a significantly different abundance in resistant cattle following early and prolonged tick exposures (compared to resistant naïve) were associated with immune response, blood coagulation, homeostasis, and wound healing. In contrast, susceptible cattle developed some of these responses only after prolonged tick exposure. Discussion: Resistant cattle were able to transmigrate immune-response related proteins towards the tick bite sites, which may prevent tick feeding. Significantly differentially abundant proteins identified in this research in resistant naïve cattle may provide a rapid and efficient protective response to tick infestation. Physical barrier (skin integrity and wound healing) mechanisms and systemic immune responses were key contributors to resistance. Immune response-related proteins such as C4, C4a, AGP and CGN1 (naïve samples), CD14, GC and AGP (post-infestation) should be further investigated as potential biomarkers for tick resistance.

Transcriptional changes in the peripheral blood leukocytes from Brangus cattle before and after tick challenge with Rhipicephalus australis.

BACKGROUND: Disease emergence and production loss caused by cattle tick infestations have focused attention on genetic selection strategies to breed beef cattle with increased tick resistance. However, the mechanisms behind host responses to tick infestation have not been fully characterised. Hence, this study examined gene expression profiles of peripheral blood leukocytes from tick-naive Brangus steers (Bos taurus x Bos indicus) at 0, 3, and 12 weeks following artificial tick challenge experiments with Rhipicephalus australis larvae. The aim of the study was to investigate the effect of tick infestation on host leukocyte response to explore genes associated with the expression of high and low host resistance to ticks. RESULTS: Animals with high (HR, n = 5) and low (LR, n = 5) host resistance were identified after repeated tick challenge. A total of 3644 unique differentially expressed genes (FDR < 0.05) were identified in the comparison of tick-exposed (both HR and LR) and tick-naive steers for the 3-week and 12-week infestation period. Enrichment analyses showed genes were involved in leukocyte chemotaxis, coagulation, and inflammatory response. The IL-17 signalling, and cytokine-cytokine interactions pathways appeared to be relevant in protection and immunopathology to tick challenge. Comparison of HR and LR phenotypes at timepoints of weeks 0, 3, and 12 showed there were 69, 8, and 4 differentially expressed genes, respectively. Most of these genes were related to immune, tissue remodelling, and angiogenesis functions, suggesting this is relevant in the development of resistance or susceptibility to tick challenge. CONCLUSIONS: This study showed the effect of tick infestation on Brangus cattle with variable phenotypes of host resistance to R. australis ticks. Steers responded to infestation by expressing leukocyte genes related to chemotaxis, cytokine secretion, and inflammatory response. The altered expression of genes from the bovine MHC complex in highly resistant animals at pre- and post- infestation stages also supports the relevance of this genomic region for disease resilience. Overall, this study offers a resource of leukocyte gene expression data on matched tick-naive and tick-infested steers relevant for the improvement of tick resistance in composite cattle.

Genome variation in tick infestation and cryptic divergence in Tunisian indigenous sheep.

BACKGROUND: Ticks are obligate haematophagous ectoparasites considered second to mosquitos as vectors and reservoirs of multiple pathogens of global concern. Individual variation in tick infestation has been reported in indigenous sheep, but its genetic control remains unknown. RESULTS: Here, we report 397 genome-wide signatures of selection overlapping 991 genes from the analysis, using ROH, LR-GWAS, XP-EHH, and FST, of 600 K SNP genotype data from 165 Tunisian sheep showing high and low levels of tick infestations and piroplasm infections. We consider 45 signatures that are detected by consensus results of at least two methods as high-confidence selection regions. These spanned 104 genes which included immune system function genes, solute carriers and chemokine receptor. One region spanned STX5, that has been associated with tick resistance in cattle, implicating it as a prime candidate in sheep. We also observed RAB6B and TF in a high confidence candidate region that has been associated with growth traits suggesting natural selection is enhancing growth and developmental stability under tick challenge. The analysis also revealed fine-scale genome structure indicative of cryptic divergence in Tunisian sheep. CONCLUSIONS: Our findings provide a genomic reference that can enhance the understanding of the genetic architecture of tick resistance and cryptic divergence in indigenous African sheep.

Proteomic profiling of plasma-derived small extracellular vesicles: a novel tool for understanding the systemic effects of tick burden in cattle.

Cattle ticks pose a significant threat to the health and profitability of cattle herds globally. The investigation of factors leading to natural tick resistance in cattle is directed toward targeted breeding strategies that may combat cattle tick infestation on the genetic level. Exosomes (EXs), small extracellular vesicles (EVs) of 50 to 150 nm diameter, are released from all cell types into biofluids such as blood plasma and milk, have been successfully used in diagnostic and prognostic studies in humans, and can provide essential information regarding the overall health state of animals. Mass spectrometry (MS) is a highly sensitive proteomics application that can be used to identify proteins in a complex mixture and is particularly useful for biomarker development. In this proof of principle study, EXs were isolated from the blood plasma of cattle (Bos taurus) with high (HTR) and low tick resistance (LTR) (n = 3/group). Cattle were classified as HTR or LTR using a tick scoring system, and EXs isolated from the cattle blood plasma using an established protocol. EXs were subjected to MS analysis in data-dependent acquisition mode and protein search performed using Protein Pilot against the B. taurus proteome. A total of 490 unique proteins were identified across all samples. Of these, proteins present in all replicates from each group were selected for further analysis (HTR = 121; LTR = 130). Gene ontology analysis was performed using PANTHER GO online software tool. Proteins unique to HTR and LTR cattle were divided by protein class, of which 50% were associated with immunity/defense in the HTR group, whereas this protein class was not detected in EXs from LTR cattle. Similarly, unique proteins in HTR cattle were associated with B-cell activation, immunoglobins, immune response, and cellular iron ion homeostasis. In LTR cattle, unique exosomal proteins were associated with actin filament binding, purine nucleotide binding, plasma membrane protein complex, and carbohydrate derivative binding. This is the first study to demonstrate that MS analysis of EXs derived from the blood plasma of HTR and LTR cattle can be successfully applied to profile the systemic effects of tick burden.

Cattle ticks are a significant burden to cattle industries globally. Current methods to treat cattle ticks are costly and inefficient in the long term. It has been noted that while some cattle may exhibit a natural resistance to ticks, others carry a heavy tick burden. The study of small extracellular vesicles, or exosomes (EXs), isolated from cattle blood plasma provides a noninvasive way of analyzing changes at the cellular level and may be of use in understanding the systemic effects of tick burden or factors leading to natural resistance. The aim of this study was to assess high (HTR) and low tick resistance (LTR) cattle identified using a tick burden scoring system by analyzing the protein content of circulating EXs via qualitative proteomics analysis. We found that a class of proteins related to defense/immunity comprised 50% of proteins unique to HTR cattle, while this protein class was not detected in proteins unique to LTR cattle. Additionally, epidermal growth factor­calcium-binding protein domains were 2-fold increased in LTR cattle compared with HTR cattle, indicating a possible mechanism for widespread metabolic change. This is the first study to employ proteomic analysis of exosomal cargo as an approach to understanding the systemic effects of tick burden in cattle.

Multiple Country and Breed Genomic Prediction of Tick Resistance in Beef Cattle.

Ticks cause substantial production losses for beef and dairy cattle. Cattle resistance to ticks is one of the most important factors affecting tick control, but largely neglected due to the challenge of phenotyping. In this study, we evaluate the pooling of tick resistance phenotyped reference populations from multi-country beef cattle breeds to assess the possibility of improving host resistance through multi-trait genomic selection. Data consisted of tick counts or scores assessing the number of female ticks at least 4.5 mm length and derived from seven populations, with breed, country, number of records and genotyped/phenotyped animals being respectively: Angus (AN), Brazil, 2,263, 921/1,156, Hereford (HH), Brazil, 6,615, 1,910/2,802, Brangus (BN), Brazil, 2,441, 851/851, Braford (BO), Brazil, 9,523, 3,062/4,095, Tropical Composite (TC), Australia, 229, 229/229, Brahman (BR), Australia, 675, 675/675, and Nguni (NG), South Africa, 490, 490/490. All populations were genotyped using medium density Illumina SNP BeadChips and imputed to a common high-density panel of 332,468 markers. The mean linkage disequilibrium (LD) between adjacent SNPs varied from 0.24 to 0.37 across populations and so was sufficient to allow genomic breeding values (GEBV) prediction. Correlations of LD phase between breeds were higher between composites and their founder breeds (0.81 to 0.95) and lower between NG and the other breeds (0.27 and 0.35). There was wide range of estimated heritability (0.05 and 0.42) and genetic correlation (-0.01 and 0.87) for tick resistance across the studied populations, with the largest genetic correlation observed between BN and BO. Predictive ability was improved under the old-young validation for three of the seven populations using a multi-trait approach compared to a single trait within-population prediction, while whole and partial data GEBV correlations increased in all cases, with relative improvements ranging from 3% for BO to 64% for TC. Moreover, the multi-trait analysis was useful to correct typical over-dispersion of the GEBV. Results from this study indicate that a joint genomic evaluation of AN, HH, BN, BO and BR can be readily implemented to improve tick resistance of these populations using selection on GEBV. For NG and TC additional phenotyping will be required to obtain accurate GEBV.

The antiviral immunity of ticks against transmitted viral pathogens.

Ticks, being obligate hematophagous arthropods, are exposed to various blood-borne pathogens, including arboviruses. Consequently, their feeding behavior can readily transmit economically important viral pathogens to humans and animals. With this tightly knit vector and pathogen interaction, the replication and transmission of tick-borne viruses (TBVs) must be highly regulated by their respective tick vectors to avoid any adverse effect on the ticks' biological development and viability. Knowledge about the tick-virus interface, although gaining relevant advances in recent years, is advancing at a slower pace than the scientific developments related to mosquito-virus interactions. The unique and complicated feeding behavior of ticks, compared to that of other blood-feeding arthropods, also limits the studies that would further elaborate the antiviral immunity of ticks against TBVs. Hence, knowledge of molecular and cellular immune mechanisms at the tick-virus interface, will further elucidate the successful viral replication of TBVs in ticks and their effective transmission to human and animal hosts.

A combined transcriptomic approach to identify candidates for an anti-tick vaccine blocking B. afzelii transmission.

Ixodes ricinus is the vector for Borrelia afzelii, the predominant cause of Lyme borreliosis in Europe, whereas Ixodes scapularis is the vector for Borrelia burgdorferi in the USA. Transcription of several I. scapularis genes changes in the presence of B. burgdorferi and contributes to successful infection. To what extend B. afzelii influences gene expression in I. ricinus salivary glands is largely unknown. Therefore, we measured expression of uninfected vs. infected tick salivary gland genes during tick feeding using Massive Analysis of cDNA Ends (MACE) and RNAseq, quantifying 26.179 unique transcripts. While tick feeding was the main differentiator, B. afzelii infection significantly affected expression of hundreds of transcripts, including 465 transcripts after 24 h of tick feeding. Validation of the top-20 B. afzelii-upregulated transcripts at 24 h of tick feeding in ten biological genetic distinct replicates showed that expression varied extensively. Three transcripts could be validated, a basic tail protein, a lipocalin and an ixodegrin, and might be involved in B. afzelii transmission. However, vaccination with recombinant forms of these proteins only marginally altered B. afzelii infection in I. ricinus-challenged mice for one of the proteins. Collectively, our data show that identification of tick salivary genes upregulated in the presence of pathogens could serve to identify potential pathogen-blocking vaccine candidates.

A genetic and immunological comparison of tick-resistance in beef cattle following artificial infestation with Rhipicephalus ticks.

Host resistance to ticks can be explored as a possible approach of combating tick infestations to complement the existing unsustainable tick control methods. Thirty-six beef cattle animals were used, consisting of Angus, Brahman and Nguni breeds, with each breed contributing 12 animals. Half of the animals per breed were artificially challenged with Rhipicephalus microplus and the other half with R. decoloratus unfed larvae per animal. Skin biopsies and blood samples were collected pre-infestation and 12 h post-infestation from the feeding sites of visibly engorging ticks. The success rate of the ticks was high and had an influence even at the early time point. Increased lymphocytes and blood urea nitrogen levels as well as decreased levels of segmented neutrophils were observed in the Angus, which were the opposite of those in the Brahman and Nguni. The increase in cholesterol, which was highest in the Angus and lowest in the Nguni, may be due to altered protein metabolism. The expression profiles of genes TRAF6, TBP, LUM and B2M were significantly different among breeds. Five genes (CCR1, TLR5, TRAF6, TBP, BDA20) had increased or constant expression post-infestation, whereas the expression of CXCL8, IL-10 and TNF-α decreased or remained the same after tick challenge. Genes that showed variation are involved in discouraging long-term supply of blood meal to the tick and those associated with immune responses. The gene LUM is a potential biomarker for tick resistance in cattle. The response to infestation by the breeds was consistent across the tick species.

SNP Diversity in CD14 Gene Promoter Suggests Adaptation Footprints in Trypanosome Tolerant N'Dama (Bos taurus) but not in Susceptible White Fulani (Bos indicus) Cattle.

Immune response to infections has been shown to be mediated by genetic diversity in pattern recognition receptors, leading to disease tolerance or susceptibility. We elucidated naturally occurring variations within the bovine CD14 gene promoter in trypanosome-tolerant (N'Dama) and susceptible (White Fulani) cattle, with genomic and computational approaches. Blood samples were collected from White Fulani and N'Dama cattle, genomic DNA extracted and the entire promoter region of the CD14 gene amplified by PCR. We sequenced this region and performed in silico computation to identify SNP variants, transcription factor binding sites, as well as micro RNAs in the region. CD14 promoter sequences were compared with the reference bovine genome from the Ensembl database to identify various SNPs. Furthermore, we validated three selected N'Dama specific SNPs using custom Taqman SNP genotyping assay for genetic diversity. In all, we identified a total of 54 and 41 SNPs at the CD14 promoter for N'Dama and White Fulani respectively, including 13 unique SNPs present in N'Dama only. The significantly higher SNP density at the CD14 gene promoter region in N'Dama may be responsible for disease tolerance, possibly an evolutionary adaptation. Our genotype analysis of the three loci selected for validation show that mutant alleles (A/A, C/C, and A/A) were adaptation profiles within disease tolerant N'Dama. A similar observation was made for our haplotype analysis revealing that haplotypes H1 (ACA) and H2 (ACG) were significant combinations within the population. The SNP effect prediction revealed 101 and 89 new transcription factor binding sites in N'Dama and White Fulani, respectively. We conclude that disease tolerant N'Dama possessing higher SNP density at the CD14 gene promoter and the preponderance of mutant alleles potentially confirms the significance of this promoter in immune response, which is lacking in susceptible White Fulani. We, therefore, recommend further in vitro and in vivo study of this observation in infected animals, as the next step for understanding genetic diversity relating to varying disease phenotypes in both breeds.

Network analysis uncovers putative genes affecting resistance to tick infestation in Braford cattle skin.

BACKGROUND: Genetic resistance in cattle is considered a suitable way to control tick burden and its consequent losses for livestock production. Exploring tick-resistant (R) and tick-susceptible (S) hosts, we investigated the genetic mechanisms underlying the variation of Braford resistance to tick infestation. Skin biopsies from four-times-artificially infested R (n = 20) and S (n = 19) hosts, obtained before the first and 24 h after the fourth tick infestation were submitted to RNA-Sequencing. Differential gene expression, functional enrichment, and network analysis were performed to identify genetic pathways and transcription factors (TFs) affecting host resistance. RESULTS: Intergroup comparisons of hosts before (Rpre vs. Spre) and after (Rpost vs. Spost) tick infestation found 51 differentially expressed genes (DEGs), of which almost all presented high variation (TopDEGs), and 38 were redundant genes. Gene expression was consistently different between R and S hosts, suggesting the existence of specific anti-tick mechanisms. In the intragroup comparisons, Rpost vs. Rpre and Spost vs. Spre, we found more than two thousand DEGs in response to tick infestation in both resistance groups. Redundant and non-redundant TopDEGs with potential anti-tick functions suggested a role in the development of different levels of resistance within the same breed. Leukocyte chemotaxis was over-represented in both hosts, whereas skin degradation and remodeling were only found in TopDEGs from R hosts. Also, these genes indicated the participation of cytokines, such as IL6 and IL22, and the activation of Wingless (WNT)-signaling pathway. A central gene of this pathway, WNT7A, was consistently modulated when hosts were compared. Moreover, the findings based on a genome-wide association study (GWAS) corroborate the prediction of the WNT-signaling pathway as a candidate mechanism of resistance. The regulation of immune response was the most relevant pathway predicted for S hosts. Members of Ap1 and NF-kB families were the most relevant TFs predicted for R and S, respectively. CONCLUSION: This work provides indications of genetic mechanisms presented by Braford cattle with different levels of resistance in response to tick infestation, contributing to the search of candidate genes for tick resistance in bovine.

Differential Tick Salivary Protein Profiles and Human Immune Responses to Lone Star Ticks (Amblyomma americanum) From the Wild vs. a Laboratory Colony.

Ticks are a growing concern to human and animal health worldwide and they are leading vectors of arthropod-borne pathogens in the United States. Ticks are pool blood feeders that can attach to the host skin for days to weeks using their saliva to counteract the host defenses. Tick saliva, as in other hematophagous arthropods, contains pharmacological and immunological active compounds, which modulate local and systemic immune responses and induce antibody production. In the present study, we explore differences in the salivary gland extract (SGE) protein content of Amblyomma americanum ticks raised in a laboratory colony (CT) vs. those collected in the field (FT). First, we measured the IgG antibody levels against SGE in healthy volunteers residing in Kansas. ELISA test showed higher IgG antibody levels when using the SGE from CT as antigen. Interestingly, antibody levels against both, CT-SGE and FT-SGE, were high in the warm months (May-June) and decreased in the cold months (September-November). Immunoblot testing revealed a set of different immunogenic bands for each group of ticks and mass spectrometry data revealed differences in at 19 proteins specifically identified in the CT-SGE group and 20 from the FT-SGE group. Our results suggest that differences in the salivary proteins between CT-SGE and FT-SGE may explain the differential immune responses observed in this study.

Genome-wide association studies for tick resistance in Bos taurus × Bos indicus crossbred cattle: A deeper look into this intricate mechanism.

Rhipicephalus (Boophilus) microplus is the main cattle ectoparasite in tropical areas. Gir × Holstein crossbred cows are well adapted to different production systems in Brazil. In this context, we performed genome-wide association study (GWAS) and post-GWAS analyses for R. microplus resistance in an experimental Gir × Holstein F2 population. Single nucleotide polymorphisms (SNP) identified in GWAS were used to build gene networks and to investigate the breed of origin for its alleles. Tick artificial infestations were performed during the dry and rainy seasons. Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA) and single-step BLUP procedure was used for GWAS. Post-GWAS analyses were performed by gene ontology terms enrichment and gene transcription factors networks, generated from enriched transcription factors, identified from the promoter sequences of selected gene sets. The genetic origin of marker alleles in the F2 population was assigned using the breed of origin of alleles approach. Heritability estimates for tick counts were 0.40 ± 0.11 in the rainy season and 0.54 ± 0.11 in the dry season. The top ten 0.5-Mbp windows with the highest percentage of genetic variance explained by SNP markers were found in chromosomes 10 and 23 for both the dry and rainy seasons. Gene network analyses allowed the identification of genes involved with biological processes relevant to immune system functions (TREM1, TREM2, and CD83). Gene-transcription factors network allowed the identification of genes involved with immune functions (MYO5A, TREML1, and PRSS16). In resistant animals, the average proportion of animals showing significant SNPs with paternal and maternal alleles originated from Gir breed was 44.8% whereas the proportion of animals with both paternal and maternal alleles originated from Holstein breed was 11.3%. Susceptible animals showing both paternal and maternal alleles originated from Holstein breed represented 44.6% on average, whereas both paternal and maternal alleles originated from Gir breed animals represented 9.3%. This study allowed us to identify candidate genes for tick resistance in Gir × Holstein crossbreds in both rainy and dry seasons. According to the origin of alleles analysis, we found that most animals classified as resistant showed 2 alleles from Gir breed, while the susceptible ones showed alleles from Holstein. Based on these results, the identified genes may be thoroughly investigated in additional experiments aiming to validate their effects on tick resistance phenotype in cattle.

Expression Patterns of Host Inflammatory Cytokine Genes during Infestation with Haemaphysalis longicornis, a Zoonotic Vector, in Blood Sucking Periods.

Tick saliva is critically important for continuous attachment to the host, blood feeding for days, and transmission of tick-borne pathogens. To characterize the patterns of inflammatory cytokine gene expression during its attachment and blood sucking time, peripheral blood samples of rabbits infested with Haemaphysalis longicornis ticks were collected at different intervals. Blood histamine concentration was evaluated as well as gene encoding IFN-γ, TNF-α, IL-2, IL-6, IL-4, and IL-10 were compared with non-infested rabbits. Blood histamine concentration of tick-infested rabbits during fast feeding time was significantly higher than that of non-infested rabbits. In both nymph and adult tick infested rabbits, expression of TNF-α and IFN-γ genes were decreased significantly (P<0.05), while expression of IL-4, IL-6, and IL-10 were increased 1.3 to 7 folds in adult infested rabbits with the exception of IL-6 that was significantly (P<0.05) decreased in nymph infested rabbits. IL-2 was not expressed in either nymph or adult infestation. H. longicornis saliva is capable of modulate host responses through a complex correlation with histamine and Th1, Th2 mediated cytokines that suppress the inflammatory responses directed toward inflammatory mediators introduced into the host during tick feeding.

Genetic parameters for tick counts across months for different tick species and anatomical locations in South African Nguni cattle.

The objective of the study was to characterise genetic parameters across months for different tick species and anatomical locations in South African Nguni cattle. Tick counts were conducted monthly, over a 2-year period, on 586 Nguni cattle under natural infestation, from four herds located in different provinces of South Africa. The counts were recorded for six species of ticks (Amblyomma hebraeum, Rhipicephalus evertsi evertsi, Rhipicephalus decoleratus and microplus (Boofilids), Rhipicephalus appendiculatus, Rhipicephalus simus and Hyalomma marginatum) attached on eight anatomical locations on the animals and were summed by species and anatomical location. Heritability estimates, phenotypic and genetic correlations were estimated on a monthly basis using mixed linear models, fitting univariate and bivariate sire models. Fixed effects considered were location, sex, year and age as a covariate. Tick counts were higher in the hot months, and A. hebraeum was the most dominant tick species. Heritability estimates for tick count varied by month and trait and ranged from 0 to 0.89. Genetic correlations were mostly positive, and low to high, with some negative correlations with high standard error. Phenotypic correlations were low to moderate. In general, high genetic correlations were observed between whole body count and the anatomical location counts, suggesting that it may not be necessary to conduct whole body counts. Counts from the belly and perineum appeared to be the most suitable surrogate traits for whole body count. These findings provide useful information for developing strategies for the practical implementation of genetic selection, as a supplement to the traditional tick control measures.

Tag SNP selection for prediction of tick resistance in Brazilian Braford and Hereford cattle breeds using Bayesian methods.

BACKGROUND: Cattle resistance to ticks is known to be under genetic control with a complex biological mechanism within and among breeds. Our aim was to identify genomic segments and tag single nucleotide polymorphisms (SNPs) associated with tick-resistance in Hereford and Braford cattle. The predictive performance of a very low-density tag SNP panel was estimated and compared with results obtained with a 50 K SNP dataset. RESULTS: BayesB (π = 0.99) was initially applied in a genome-wide association study (GWAS) for this complex trait by using deregressed estimated breeding values for tick counts and 41,045 SNP genotypes from 3455 animals raised in southern Brazil. To estimate the combined effect of a genomic region that is potentially associated with quantitative trait loci (QTL), 2519 non-overlapping 1-Mb windows that varied in SNP number were defined, with the top 48 windows including 914 SNPs and explaining more than 20% of the estimated genetic variance for tick resistance. Subsequently, the most informative SNPs were selected based on Bayesian parameters (model frequency and t-like statistics), linkage disequilibrium and minor allele frequency to propose a very low-density 58-SNP panel. Some of these tag SNPs mapped close to or within genes and pseudogenes that are functionally related to tick resistance. Prediction ability of this SNP panel was investigated by cross-validation using K-means and random clustering and a BayesA model to predict direct genomic values. Accuracies from these cross-validations were 0.27 ± 0.09 and 0.30 ± 0.09 for the K-means and random clustering groups, respectively, compared to respective values of 0.37 ± 0.08 and 0.43 ± 0.08 when using all 41,045 SNPs and BayesB with π = 0.99, or of 0.28 ± 0.07 and 0.40 ± 0.08 with π = 0.999. CONCLUSIONS: Bayesian GWAS model parameters can be used to select tag SNPs for a very low-density panel, which will include SNPs that are potentially linked to functional genes. It can be useful for cost-effective genomic selection tools, when one or a few key complex traits are of interest.

In silico target network analysis of de novo-discovered, tick saliva-specific microRNAs reveals important combinatorial effects in their interference with vertebrate host physiology.

The hard tick Ixodes ricinus is an important disease vector whose salivary secretions mediate blood-feeding success on vertebrate hosts, including humans. Here we describe the expression profiles and downstream analysis of de novo-discovered microRNAs (miRNAs) expressed in I. ricinus salivary glands and saliva. Eleven tick-derived libraries were sequenced to produce 67,375,557 Illumina reads. De novo prediction yielded 67 bona fide miRNAs out of which 35 are currently not present in miRBase. We report for the first time the presence of microRNAs in tick saliva, obtaining furthermore molecular indicators that those might be of exosomal origin. Ten out of these microRNAs are at least 100 times more represented in saliva. For the four most expressed microRNAs from this subset, we analyzed their combinatorial effects upon their host transcriptome using a novel in silico target network approach. We show that only the inclusion of combinatorial effects reveals the functions in important pathways related to inflammation and pain sensing. A control set of highly abundant microRNAs in both saliva and salivary glands indicates no significant pathways and a far lower number of shared target genes. Therefore, the analysis of miRNAs from pure tick saliva strongly supports the hypothesis that tick saliva miRNAs can modulate vertebrate host homeostasis and represents the first direct evidence of tick miRNA-mediated regulation of vertebrate host gene expression at the tick-host interface. As such, the herein described miRNAs may support future drug discovery and development projects that will also experimentally question their predicted molecular targets in the vertebrate host.

Design of selection schemes to include tick resistance in the breeding goal for Hereford and Braford cattle.

Ticks are one of the main causes of losses in cattle, causing economic impact by reducing productivity and fertility and by transmission of diseases. The objective of this study was to analyze the genetic gains obtained through different strategies to include traditional (EBV) or genomic EBV (GEBV) for tick count (TC) in selection indexes for Hereford and Braford cattle. Besides TC, we also considered traits currently included in the Delta G Breeding Program Index (DGI): preweaning gain, weaning conformation, weaning precocity, weaning muscling, postweaning gain, yearling conformation, yearling precocity, yearling muscling, and scrotal circumference. Genetic gain per generation (ΔG) was evaluated using the current DGI and including TC in 8 alternative scenarios with TC relative weightings of 10, 50, or 100% and using phenotype or GEBV. Genomic EBV accuracy () ranged between 0.1 and 0.9. As expected, increasing increases the accuracy of the index () for all scenarios in which GEBV were considered. As the relative weight of TC was increased to 50%, greater ΔG differences in relation to the baseline DGI ($53.03) scenario were observed when the GEBV information was included with equal to or greater than 0.7 only for TC (ΔG between $61.06 and $74.26) or equal to or greater than 0.5 for all traits (ΔG between $56.03 and $83.36). To achieve these accuracies for traits with low heritability, a large calibration data set would be required. Focusing only on TC, the availability of genomic information would be desirable to avoid the need to count ticks and the exposure of animals to parasitism risks. However, for = 0.7, the respective numbers for Hereford and Braford would be 4,703 and 6,522 animals. As expected, when comparing the relative index weights of 10, 50, and 100% for TC, the highest response to selection per generation (RS) for TC was in the scenario was with 100% relative weight and GEBV for this trait (SR = -0.09 SD with = 0.9). This would be the recommended scenario to form tick-resistant lines in Hereford and Braford cattle. However, with 50% relative weight for TC, including GEBV information for TC only or for all traits in index ( = 0.9), it should yield 93 or 84% of RS, respectively, compared to that obtained with full emphasis on TC (100% relative weight) and GEBV information. This indicates that in the presence of highly accurate GEBV, despite slightly slower gain for TC, indexes with 50% relative weight for TC are interesting alternatives to jointly improve tick resistance and other relevant traits.

Immune recognition of salivary proteins from the cattle tick Rhipicephalus microplus differs according to the genotype of the bovine host.

BACKGROUND: Males of the cattle tick Rhipicephalus microplus produce salivary immunoglobulin-binding proteins and allotypic variations in IgG are associated with tick loads in bovines. These findings indicate that antibody responses may be essential to control tick infestations. Infestation loads with cattle ticks are heritable: some breeds carry high loads of reproductively successful ticks, in others, few ticks feed and they reproduce inefficiently. Different patterns of humoral immunity against tick salivary proteins may explain these phenotypes. METHODS: We describe the profiles of humoral responses against tick salivary proteins elicited during repeated artificial infestations of bovines of a tick-resistant (Nelore) and a tick-susceptible (Holstein) breed. We measured serum levels of total IgG1, IgG2 and IgE immunoglobulins and of IgG1 and IgG2 antibodies specific for tick salivary proteins. With liquid chromatography followed by mass spectrometry we identified tick salivary proteins that were differentially recognized by serum antibodies from tick-resistant and tick-susceptible bovines in immunoblots of tick salivary proteins separated by two-dimensional electrophoresis. RESULTS: Baseline levels of total IgG1 and IgG2 were significantly higher in tick-susceptible Holsteins compared with resistant Nelores. Significant increases in levels of total IgG1, but not of IgG2 accompanied successive infestations in both breeds. Resistant Nelores presented with significantly higher levels of salivary-specific antibodies before and at the first challenge with tick larvae; however, by the third challenge, tick-susceptible Holsteins presented with significantly higher levels of IgG1 and IgG2 tick salivary protein-specific antibodies. Importantly, sera from tick-resistant Nelores reacted with 39 tick salivary proteins in immunoblots of salivary proteins separated in two dimensions by electrophoresis versus only 21 spots reacting with sera from tick-susceptible Holsteins. CONCLUSIONS: Levels of tick saliva-specific antibodies were not directly correlated with infestation phenotypes. However, in spite of receiving apparently lower amounts of tick saliva, tick-resistant bovines recognized more tick salivary proteins. These reactive salivary proteins are putatively involved in several functions of parasitism and blood-feeding. Our results indicate that neutralization by host antibodies of tick salivary proteins involved in parasitism is essential to control tick infestations.

Immune and biochemical responses in skin differ between bovine hosts genetically susceptible and resistant to the cattle tick Rhipicephalus microplus.

BACKGROUND: Ticks attach to and penetrate their hosts' skin and inactivate multiple components of host responses in order to acquire a blood meal. Infestation loads with the cattle tick, Rhipicephalus microplus, are heritable: some breeds carry high loads of reproductively successful ticks, whereas in others, few ticks feed and reproduce efficiently. METHODS: In order to elucidate the mechanisms that result in the different outcomes of infestations with cattle ticks, we examined global gene expression and inflammation induced by tick bites in skins from one resistant and one susceptible breed of cattle that underwent primary infestations with larvae and nymphs of R. microplus. We also examined the expression profiles of genes encoding secreted tick proteins that mediate parasitism in larvae and nymphs feeding on these breeds. RESULTS: Functional analyses of differentially expressed genes in the skin suggest that allergic contact-like dermatitis develops with ensuing production of IL-6, CXCL-8 and CCL-2 and is sustained by HMGB1, ISG15 and PKR, leading to expression of pro-inflammatory chemokines and cytokines that recruit granulocytes and T lymphocytes. Importantly, this response is delayed in susceptible hosts. Histopathological analyses of infested skins showed inflammatory reactions surrounding tick cement cones that enable attachment in both breeds, but in genetically tick-resistant bovines they destabilized the cone. The transcription data provided insights into tick-mediated activation of basophils, which have previously been shown to be a key to host resistance in model systems. Skin from tick-susceptible bovines expressed more transcripts encoding enzymes that detoxify tissues. Interestingly, these enzymes also produce volatile odoriferous compounds and, accordingly, skin rubbings from tick-susceptible bovines attracted significantly more tick larvae than rubbings from resistant hosts. Moreover, transcripts encoding secreted modulatory molecules by the tick were significantly more abundant in larval and in nymphal salivary glands from ticks feeding on susceptible bovines. CONCLUSIONS: Compared with tick-susceptible hosts, genes encoding enzymes producing volatile compounds exhibit significantly lower expression in resistant hosts, which may render them less attractive to larvae; resistant hosts expose ticks to an earlier inflammatory response, which in ticks is associated with significantly lower expression of genes encoding salivary proteins that suppress host immunity, inflammation and coagulation.